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Abstract 

The two-dimensional Ising model is the simplest model of statistical mechanics exhibiting a 
second order phase transition. While in absence of magnetic field it is known to be solvable 
on the lattice since Onsager's work of the forties, exact results for the magnetic case have 
been missing until the late eighties, when A. Zamolodchikov solved the model in a field at 
the critical temperature, directly in the scaling limit, within the framework of integrable 
quantum field theory. In this article we review this field theoretical approach to the Ising 
universality class, with particular attention to the results obtained starting from Zamolod- 
chikov's scattering solution and to their comparison with the numerical estimates on the 
lattice. The topics discussed include scattering theory, form factors, correlation functions, 
universal amplitude ratios and perturbations around integrable directions. Although we re- 
strict our discussion to the Ising model, the emphasis is on the general methods of integrable 
quantum field theory which can be used in the study of all universality classes of critical 
behaviour in two dimensions. 



1 Introduction 



A lattice system close to a second order phase transition point exhibits a number of features 
which do not depend on the specific microscopic realisation and coincide with those of all systems 
sharing the same essential symmetry properties in the given spatial dimension. In principle, the 
continous field theoretical description of the scaling region provides the most natural theoretical 
framework for the quantitative study of these universal features; in practice, however, the need 
of non-perturbative methods seriously complicates the task. 

The Ising model is the fundamental model in the theory of critical phenomena. Its theo- 
retical importance became evident in 1944, when Onsager was able to compute the free energy 
on the square lattice in absence of magnetic field, providing in this way the first exact description 
of a second order phase transition 2 . It had to become clear later that the Ising model corre- 
sponds to the simplest universality class of critical behaviour. Since Onsager's work, the Ising 
model has been an essential indicator of the progress in the analytic study of critical phenomena. 
While it remains unsolved in three dimensions, the lattice studies gave additional exact results 
in the two dimensional case. In 1952 Yang published the first derivation of the formula for the 
spontaneous magnetisation that Onsager had presented three years before [3]. Scaling theory 
would later show that this result amounted to completing the list of critical exponents of the 
Ising universality class. Further important progress was made on the determination of correla- 
tion functions .4\ It was found, in particular, that the spin-spin correlator can be expressed, in 
the scaling limit, through the solution of a differential equation of Painleve type 0. Up to some 
generalisations, this remains the only non-trivial correlation function of quantum field theory to 
be exactly known. 

These results refer to zero magnetic field, the Ising model in a field having never been solved 
on the lattice. For this reason, A. Zamolodchikov's solution of the model with magnetic field, 
at the critical temperature, directly in the scaling limit, came as a major surprise for many 
Even more so since this solution consisted of a long list of scattering amplitudes for eight 
different species of relativistic particles. The most striking aspect of Zamolodchikov's work, 
however, was that, apart from its consequences for the Ising model, it actually implied new 
exact results for all the universality classes of critical behaviour in two dimensions. In fact, if 
conformal field theory had given a complete description of critical points 0, now it was shown 
that some renormalisation group trajectories flowing out of each critical point correspond to 
exactly solvable (integrable) quantum field theories. These describe particular directions in 
the scaling region of statistical models which, in general, are not solvable as long as the non- 
universal lattice details are not eliminated through the scaling limit. Specific realisations of a 
given universality class, however, can be solvable already on the lattice. In particular, a solvable 
lattice model yielding the same scaling limit of the Ising model in a magnetic field at critical 
temperature was found in 8 . 

Dealing directly with the scaling limit, integrable quantum field theory is the most effective 
tool for extracting exact information about universality classes. In this context, 'integrable' 
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means that the relativistic scattering theory associated to the quantum field theory can be 
determined exactly. Once this has been done, the next task is that of bridging the gap between 
the scattering solution and the quantities of more direct interest for statistical mechanics. In 
this article we review the results obtained through this approach for the universality class of 
the Ising model in a magnetic field, on the infinite plane. Although we refrain from making 
reference to other models, we stress the general reasons which allow to obtain similar results for 
the other universality classes in two dimensions. Since we always work in the continuum limit, 
we take care of comparing the field theoretical predictions for the universal quantities with the 
available lattice estimates. 

The article is organised as follows. In the next section we recall the definition of the Ising 
model on the lattice before turning to the field theoretical description of the critical point 
and the scaling region around it. In section 3 we review the origin of integrable quantum 
field theories, their solution in the scattering framework and the form factor approach to the 
computation of correlation functions. All this is illustrated in practice through the application 
to the integrable directions of the Ising field theory in section 4, where the known results for 
the purely thermal case and the recent advances for the magnetic case are presented within the 
same framework. In particular, we discuss the determination of form factors in the magnetic case 
Elj • Section 5 illustrates how the recent field theoretical results reflect onto the traditional 
way of characterising critical behaviour and allow to complete the list of canonical amplitude 
ratios for the Ising universality class. Finally, in section 6, we briefly discuss how to exploit the 
integrable directions for a more general investigation of the scaling region and, in particular, 
review few basic results on the evolution of the particle spectrum in the Ising field theory. 



2 Ising field theory in two dimensions 

The Ising model is defined on the lattice by the reduced Hamiltonian 

tf = ~5><°-i-ir£< 7 *> (2.i) 

(iJ) i 

where Oi = ±1 is the spin variable at the i-th site and the first sum is taken over nearest neigh- 
bours; the couplings T and H are referred to as temperature and magnetic field, respectively. 
The expectation value of any lattice variable O is given by 

(0) = W° e ~ E i (2-2) 

where 

Z=Y,e~ E (2.3) 

R} 

is the partition function. We will always refer to the two-dimensional ferromagnetic (T > 0) 
case in the following. 
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Consider the case H = 0. The Hamiltonian is then invariant under the change of sign of all 
spins. This spin reversal symmetry is broken spontaneously when T is smaller than a critical 
value T c for which a second order phase transition takes place [2]. The critical point divides the 
temperature axis into a high-temperature, disordered phase, and a low-temperature phase where 
a spontaneous magnetisation exists. The two phases are related by a duality transformation 11. 

No global symmetry is left in the model when H ^ 0, and that located at (T, H) = (T c , 0) 
in the T-H plane is the only critical point in the model. Since changing the sign of H simply 
amounts to a global spin reversal transformation, it is sufficient to refer to the case H > 0. 

The correlation length £ diverges at a second order phase transition point and remains much 
larger than the lattice spacing in a neighbourhood of this point in coupling space. In the scaling 
region in which T — ► T c , H — > 0, the system can effectively be considered as translationally and 
rotationally invariant, and quantum field theory provides a continous description suitable for 
the investigation of the universal properties (see e.g. |12j). 

In particular, the behaviour of the critical point on scales much larger than the lattice spacing 
is described by a massless (mass ~ l/£) field theory invariant under scale transformations and, 
actually, under the larger group of conformal transformations. In two dimensions conformal 
symmetry is infinite dimensional and, for this reason, conformal field theories are exactly solved 
0. As it should be, the Ising critical point is described by the simplest (i.e. with the smallest 
operator content) conformal field theory satisfying the requirement of reflection positivity |13j . 
This theory contains three fundamental (primary) operators which are invariant (scalar) under 
rotations, namely the identity /, the spin a(x) and the energy e{x) (x = (xi, x^) denotes a point 
on the plane). The spin and energy operators are the continous version of the lattice variables 
0"i and J2j cr i cr j (j nearest neighbour of i), respectively. Each primary operator possesses an 
infinite number of 'descendents', the simplest example being provided by the derivatives of the 
primaries. A primary and its descendents share the same internal symmetry properties and form 
a 'conformal family'. 

In a conformal field theory the product of two operators and $2 (to be thought inside 
a correlation function) can be expanded over a complete basis made of an infinite number of 
operators in the form 

*i(x)$ 2 (0) = E C& a z A ^ ^ A *i - A * 2 " A *x ~ A * 2 A k (0) . (2.4) 

k 

Here z = x\ + 1x2 and z = x\ — ix2 are complex coordinates, the C^^ 2 's are called structure 
constants and A$ and A<j> are the conformal dimensions of an operator $(x). The scaling 
dimension 

A$ = A$ + A$ (2.5) 

and the euclidean spin 

Scj> = A<j> — A$ (2.6) 
determine the behaviour of the operator under scale transformations and rotations, respectively. 
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The operators with a definite scaling dimension are called scaling operators. A scalar operator 
has s<j> = 0. 

While the scaling dimension of the identity operator vanishes, the two non-trivial primary 
operators in Ising field theory have dimensions X a = 1/8 and X e = 1. These two numbers 
determine all the critical exponents of the Ising universality class (see section 5). The conformal 
dimensions of the descendent operators differ by positive integers from those of the corresponding 
primary. 

The structure of the operator product expansion in the Ising conformal theory can be sym- 
bolically expressed as 

(TX(T~[/] + [e] 

a x e ~ [<r] (2.7) 

£ X £ ~ [I] , 

where the square brackets indicate the appearence of a whole conformal family on the r.h.s. 
Clearly, this structure is compatible with the fact that the conformal families [I] and [e] are 
even under spin reversal while [a] is odd. 

The Ising critical point is left invariant by the duality transformation which exchanges the 
high- and low-temperature phases at H = 0. While the energy operator e(x), which drives the 
model away from criticality along the thermal axis, changes sign under duality (this is why [e] 
does not appear in the last of (|2.7j0 . such a transformation maps c(x) onto a disorder operator 
/i(x) |14| . The operators a and \i have the same scaling dimension 1/8, but are mutually non- 
local, in the following sense. Two operators <I>i and $2 are said to be mutually local if their 
product is unchanged when one of them is taken once around the other on the plane (i.e. under 
the analytic continuation z — > e 2t7T z, z — > e~ 2t7T z in (|2.4j0 . Of course this amounts to a statement 
about the single-valuedness of correlation functions involving the two operators. An operator is 
said to be local if it is local with respect to itself. It can be shown that the local operators (the 
only ones of interest for us) are those with integer or half-integer euclidean spin. The mildest 
type of mutual non-locality (called semi-locality) corresponds to the case 

(• • • ^(e 2m z, e~ 24 ^)d> 2 (0) • • •) = h^ 2 (- ■ • z)<D 2 (0) • • •) , (2.8) 

where phase called semi-locality factor. 

The leading term in the x — > expansion of the product a x fi is |14j 

ct(xHO) ~ \x\-V 4 [s/zil>(0) + v^^(0)] + . . . , (2.9) 

where ip and ip have conformal dimensions (A, A) = (1/2, 0) and (0, 1/2), respectively (the terms 
omitted in the r.h.s. contain descendents of if) and V>). This means that taking a around \x once 
produces a minus sign, so that the two operators are semi-local with = — 1. Similarly, ip and 
i[) are semi- local with respect to a and /i with the same semi- locality factor —1. 

Summarising, the three conformal families [/], [a] and [e] we originally considered form a 
'local section' of the Ising conformal field theory, namely a maximal set of fields all mutually 
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Figure 1: Coupling space of the Ising field theory ()2.12|) . The oriented lines indicate some 
renormalisation group trajectories flowing out of the critical point at the origin. The integrable 
directions coincide with the principal axes. In our conventions r > corresponds to T > T c . 

local and closed under operator product expansion. Duality leads to consider a second local 
section differing from the first one by the substitution of [a] with \p]. A last local section is 
given by [I], [ip], [i/j] and [e]. All the conformal families together form the full space of local 
operators of the Ising conformal theory. 

The operators ip and tp are identified by their conformal dimensions as the two components 
of a neutral (Majorana) fermion. This is why the Ising conformal point is described by the free 1 
Hamiltonian (or euclidean action) 



where d = d z = {d\ — id2)/2 and d = dz = {d\ + id2)/2. The conformal dimensions (1/2, 1/2) 
of the energy operator show that it is bilinear in the fermion components: 



a and on the contrary, are semi-local with respect to ip and if) and give rise to the non-trivial 
sector of the theory ()2.1(Jj) . 

The field theory describing the scaling region around the critical point is obtained by adding 
to the conformal action ()2.1U|) the contributions of the operators conjugated to the termperature 
and the magnetic field, namely the energy and the spin operator, respectively. This leads to the 
Ising field theory 




(2.10) 



(2.11) 




(2.12) 



where 



r ~ M 2 ~ x - = M 



1 No interaction term preserving scale invariance can be formed. 
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h ~ M 2 ~ x ° = M 15/8 



are dimensional couplings measuring the deviation from critical temperature and the magnetic 
field, respectively. Here, M ~ l/£ denotes the mass scale associated to the breaking of scale 
invariance away from criticality. It is worth stressing that the action IJ2.12J) is uniquely selected 
as the scaling limit of (|2.1j) by the fact that a and e are the only (non-constant) scalar rele- 
vant 2 operators in the local section of the operator space containing a. Additional terms in 
()2.12j) play a role only when trying to account for subleading terms in the expansion of lattice 
observables around the critical point (corrections to scaling). The field theory (|2.12j) describes 
a one-parameter family of renormalisation group trajectories flowing out of the critical point at 
r = h = (Figure 1) and labelled by the dimensionless quantity 

71 = WnE ■ ( 2 - 13 ) 

The Ising model with if = is solvable on the lattice and then must be solvable in the 
scaling limit. In fact, (|2.1()|) and (|2.11|) imply that for h = the action (|2.12|) describes a free 
massive fermion, the mass being proportional to |r|. For H ^ the Ising model has never been 
solved on the lattice. A. Zamolodchikov showed that it is solvable directly in the scaling limit 
if T = T c 6 . This is a consequence of the fact that ()2.12j) with r = is an integrable quantum 
field theory. We summarise in the next section some generalities about integrable field theories 
before turning to the study of the integrable directions in the scaling Ising model (of course the 
free theory resulting from (|2.12j) when h = is a particularly simple case of integrability) . 



3 Integrable quantum field theories 

3.1 Conserved currents 

The notion of integrability is generally associated to the presence of an infinite number of 
conserved quantities. In two-dimensional quantum field theory a conservation law takes the 
form 

Bt s+1 = ae s _i , (3.i) 

where T s+ \ and O s _i are local operators (currents) with spin s + 1 and s — 1, respectively. 
Any quantum field theory possesses the conservation law (|3.1j) with s = 1, T2 and Go being 
components of the energy-momentum tensor. Are there theories allowing for additional, non- 
trivial conservation laws? The answer is obviously affirmative for the conformal theories: any 
descendent T s of the identity with conformal dimensions (s, 0) is a local operator satisfying (13.111 
with zero on the r.h.s. Of course, this is a direct consequence of the infinite dimensional character 
of conformal symmetry in two dimensions. It is then natural to ask whether any conservation 
law (other than that of energy and momentum) can survive when a non-scale-invariant theory is 

2 In the renormalisation group language, an operator is relevant if its scaling dimension is smaller than the 
space dimensionality (2 in our case). 
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obtained as a perturbation of a conformal action Ac ft by a relevant operator <£, namely when 
considering the action 



where R]!_i are operators with conformal dimensions (s — n(l — A), 1 — n(l — A)) (the perturbing 
operator <3? has dimensions (A, A), A < 1). It follows that in any theory with a spectrum of 
conformal dimensions bounded from below the r.h.s. of 1)3.3(1 can only contain a finite number 

(n) 

of terms. Moreover, the operators R s _i with n > 1 can be accomodated within the operator 
space of a theory with a discrete spectrum of conformal dimensions only if special relations 
between the dimensions are fulfilled. Hence, in a generic case, the r.h.s. of 1)3.3(1 contains only 
the operator R^li, which is a descendent of the perturbing operator <]?. Thus the issue of 
conservation away from criticality is reduced to establishing under which conditions (if any) the 
operator gR^\ can be written in a total derivative form <90 s _2- The complete characterisation 
of the operator space provided by conformal field theory 3 allows the identification of sufficient 
conditions through a so-called 'counting argument' This exploits notions of conformal field 
theory which are not required in the remainder of this article, and we prefer to directly state the 
result of the analysis, which is remarkable: several non-trivial conservation laws (expected to be 
the first representatives of infinite series) can be found for a number of different perturbations 
of essentially all the known conformal points in two-dimensions; each integrable direction in 
coupling space is characterised by a specific set of values of the spin s for which a conservation 
law of the form ()3.1() is present. 

For the Ising field theory 1)2.12(1 . in particular, the counting argument implies integrability 
when r = and, of course, when h = 0. The theory is not integrable when both couplings are 
different from zero |151 116j . 

3.2 Scattering theory 

The field theory obtained by perturbing a conformal theory by one or more relevant operators 
normally develops a finite correlation length and admits a description in terms of massive par- 
ticles. These particles propagate in a space with one spatial and one time dimension, related 
by analytic continuation to imaginary time to the euclidean plane we have in mind for the 
applications to equilibrium statistical mechanics. 

Integrability induces major simplifications in the scattering of relativistic particles (see 
and references therein). Call 




(3.2) 



In this case the original (conformal) conservation laws get modified into [5] 



BT s =gR^ 1 + ...+g n R ( ;j 1 + ... 



(3.3) 




(3.4) 



3 It is generally assumed that the operator spaces at the conformal point and in the perturbed theory (IM.'it are 
isomorphic. 
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(a) (b) 

Figure 2: Space-time diagrams for a three-particle scattering process. A genuine three-body 
collision (a) or a sequence of widely separated two-body collisions (b) can be obtained by a 
suitable choice of initial conditions on the wave packets. 

the spin s > conserved quantities, P\ being the sum of energy and momentum. Also, denote 
by A a {pV") a particle of type a and energy-momentum p^ = (jp,p l ) satisfying the mass shell 
condition 

P% = (P ) 2 -(P X ) 2 =^- (3-5) 

A conserved quantity acts as 

Ps\A a M)---AaM)hn(out) = (f> S 0fc (Pfc)) \A ai (Pl)---AaM)in(out) (3.6) 

on the initial (final) state of a scattering process containing n widely separated particles. The 
behaviour under euclidean rotations fixes the form of the one-particle eigenvalue to be 

<4(^) = <p s , (3-7) 

with p = p° + p 1 and a constant. Then conservation means that 

n m 

E^ fe (p£) = (3.8) 
k=i j=i 

in a scattering process with n particles A ak (p^,) in the initial state and m particles Af,. (q^) 
in the final state. Since this is a system of an infinite number of equations (one for each 
conserved quantity P s ) for a finite number of unknowns (the energies and momenta of the 
outgoing particles), one concludes that in a scattering process of an integrable quantum field 
theory 

i) the final set of energies and momenta coincides with the initial one. 

This is a remarkable simplification of the relativistic scattering problem. There is, however, 
an additional important result whose origin can be understood through the following euristic 
argument 4 ^Hj. The conserved quantities P s can be seen as generators of space-time displace- 
ments on wave packets. While P\ simply translates the trajectories of all particles by the same 
4 See |19j for a discussion in the axiomatic framework. 
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Figure 3: Space-time diagram associated to the scattering amplitude S a b{6\ — 82); time runs 
upwards. 

amount in space-time, for s > 1 the amount of the shift depends on the momentum of each 
particle. Consider now a state containing three particles (wave packets) with different momenta 
and initial conditions chosen in such a way that they give rise to a genuine three-body interaction 
(Figure 2a). We can perform on this state a suitably chosen transformation generated by the 
P s with s > 1 which shifts each trajectory by a different amount and resolves the three-body 
interaction into a sequence of three two-body collisions widely separated in space-time (Figure 
2b). Since the P s are conserved, the transformation commutes with the time evolution and 
then the scattering amplitudes for the two processes coincide. The general conclusion is that in 
integrable quantum field theory 

ii) any n-particle scattering amplitude factorises into the product of n(n — l)/2 two-particle 
amplitudes. 

We discussed the elasticity (i.e. absence of particle production) and factorisation of the 
scattering assuming the presence of an infinite number of conserved quantities. It has been 
shown, however, that the existence of one such a quantity besides energy-momentum is sufficient 
to arrive at the same conclusions |2fl| . 

As a consequence of factorisation, the determination of the S'-matrix (i.e. the collection of 
all scattering amplitudes) in integrable quantum field theories reduces to that of the two-particle 
amplitudes. For the sake of simplicity we restrict our discussion to the case of neutral particles, 
which is enough for dealing with the Ising model. Then the two-particle amplitudes can be 
defined through the relation 



for a particle of mass m a . Particles are ordered with rapidities decreasing (increasing) from 
left to right in initial (final) states . The dependence of the amplitudes on rapidity differences 
is a consequence of Lorentz invariance. A pictorial representation of a two-particle scattering 




(3.9) 



where we introduced the rapidity variables parameterising on-shell momenta as 




(3.10) 
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Figure 4: Simple pole diagram associated to equation (j3,13j) . 



amplitude is shown in Figure 3. Double application of Q3.9|) together with S a b(9) = S ba (9) yield 
the unitarity relation 

Sab(9)Sab(-e) = 1 . (3.11) 

The crossing symmetry relation 

Sab(e) = S ab (i7T - 9) (3.12) 

is a general property of relativistic scattering |21| . 

The scattering amplitudes S a b(9) are meromorphic functions of the rapidity difference. A 
simple pole with residue 

TO 2 



S, 



i u: 



ah 



and u c ah G (0, ir) corresponds to a bound state with mass square 

m 2 c = m 2 a + ml + 2m a m\, cos u° ab 



(3.13) 



(3.14) 



propagating in the A a A b scattering channel. A three-particle coupling T^ b is associated to each 
vertex of the corresponding diagram in Figure 4. It follows from (|3.14j) that 



u nh = TV - u 



ab 



(3.15) 



is the angle opposite to m c in a triangle with sides of length m a , m b , m c . This implies the 
relation 

u c ab + u a bc + u b ca = 2tt . (3.16) 



As a consequence of (|3.12|) . a crossed channel pole with negative residue at 9 = 
along with each pole ()3.13j) . 

The general meromorphic solutions of (|3.11|) and (|3.12|) are of the form [22 

S ab (9) = ± [] ta(0), 



- U ab) §° eS 



(3.17) 
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Figure 5: Pictorial representation of the bootstrap equation (|3.21j) . 

with 

tanh h(6 + ma) 
ta(0) = — -fe . • (3.18) 
tann ±{t) — nra) 

Real values of a between and 1 correspond to the bound state poles ()3.13j) . Important con- 
straints on the bound state structure come from the fact that an infinite number of quantities 
besides energy and momentum has to be conserved at each three-particle vertex (Hj. In the 
vicinity of a bound state pole we can write 

+ «& - e)A b {6 - iu a bc + e)) ~ - \A C (9)) . (3.19) 

Applying P s to both sides and equating the eigenvalues gives 

K>y s# ™ + K b s m s b e- isa ^ = K c s m s c . (3.20) 

Integrability in presence of bound states also provides functional relations between different 
amplitudes. Indeed, we can consider the scattering of a particle 4 with a resonant pair 44 
and exploit the freedom of translating the particle trajectories to obtain (Figure 5) 

S dc (e) = S da (6 - iu b ac )S db (0 + iu a bc ) . (3.21) 

These equations are known as 'bootstrap' equations for the following reason. Suppose that the 
use of (|3.2fl|) plus other considerations ends up in an educated guess for the scattering amplitudes 
of the 'elementary' (lightest) particles of the theory. Then one can use (TOTT) with 4,4, 4z 
'elementary' particles and 4 ° ne of their bound states to determine the amplitude Sd c - Each 
amplitude computed in this way may present poles corresponding to new particles, in which case 
the use of (|3.21|) is iterated. Once this bootstrap procedure consistently reaches a point where 
no new particles are generated, one is left with the exact solution of an integrable quantum field 
theory (see ^3] for a series of examples). 

3.3 Correlation functions 

Having the S'-matrix is not enough for statistical mechanical applications; correlation functions 
are needed. These can be expressed in the form of spectral sums over complete sets of particle 
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Figure 6: Pictorial representation of an n-particle form factor. 



states. For example, the (euclidean) correlator of two scalar operators reads 5 



<*i(x)* 2 (0)> = E 



" (27r) n Je 1> ... > e 



where 



d6i ■ ■ ■ d6 n F a ^_ Mn (61,..., r , 



E n = Y^ m ak cosh 8 k , 
k=i 



F ai 2 ..a n (0i, • • • A) 



and the matrix elements (Figure 6) 

Fa^aJOl, • • • A) = <0|*(0)|4, a (Ox) • • • A an (9 n )} 



-\x\E n 
(3.22) 

(3.23) 
(3.24) 



are called form factors (|0) is the vacuum state). 

The form factors can be computed exactly in integrable quantum field theory. They are 
subject to a number of equations [2312] which for the case of neutral particles we are discussing 
read 



(*i + A, ... A + A) = e sA F^ an (9 u . .. A) 

F a 1 ...a t a t+1 ...a n (6l, • • • A #1+1) • • • A) = 



(3.25) 



5, 



-l)-^fi...a i+ ia ! ...a n (^l) ■ ■ ■ >#i+lA • • • A) (3.26) 

(3.27) 



R es ^-0 6 =m^FX ll ...a n (#aAA, ■ ■ ■ > #71 ) = ^fc -^cai ...a n (#c A , • • • , #n) 



+ 2iTT,9 2 , ■ ■ ■ A) = ^a/al.a^!^, • • • A A) 



(3.28) 



Res0/=e +i7r i^ ai an (9',0,0i, . . . A 



««6 ^-^n^(%- )J ^...a^l,-" A) • (3.29) 

The first equation expresses Lorentz covariance for a spin s operator, while the second 
is a consequence of (|3.9|) and factorisation of multi-particle scattering. The third equation 
5 Particle states are normalised through the condition (A a (6i)\Ab(03)) = 27v5 a tS(6i — 62). 
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says that the form factors inherit the direct channel bound state poles of the S'-matrix [9 C = 

(ut c e a + u h ca e b )/ui b ). 

Within the framework of relativistic scattering theory, a crossing process corresponds to 
trading a particle of energy-momentum p^ in the initial (final) state with an antiparticle of 
energy-momentum — p^ in the final (initial) state [2I| • The inversion of energy-momentum cor- 
responds to an in shift in the rapitidy parameterisation. Hence, equation (|3.28[) shows that the 
(kinematically immaterial) double crossing of a particle corresponds to a reordering of rapidities. 
If the operator <j) a which creates the particle 6 A a is non-local (semi-local, for simplicity) with 
respect to <J>, the double crossing of this particle also produces a phase factor l$ j( p a ^ 1 (recall 
the definition (J2SJ)) |25|. 

Finally, equation (|3.29|) expresses the fact that two identical particles with opposite energy- 
momentum can annihilate. This gives rise to a pole whose residue reflets the fact that the 
annihilation of adjacent particles can take place either directly or through the analytic contin- 
uation (JH2EJ). 

The form factor equations listed above are satisfied by all the operators of the theory and 
then admit an infinite number of solutions. Identifying the solution corresponding to a given 
operator is a central problem in this approach. The following argument turns out to be quite 
helpful in this respect. Consider a scalar operator <3? with scaling dimension X$ in a massive 
theory. Since 

($(x)*(0))~|^p-, |a;|-0 (3.30) 

we know that 7 

M p = J d 2 x\x\ p (<S>(x)<Z>(0)) c (3.31) 

is finite only if 

p + 2 > 2Xq, . (3.32) 

Once we use ()3.22|) to expand the correlator in (|3.31|) and perform the spatial intergration, we 
are left with a series of n-fold integrals over rapidities with positive integrands (moduli squared 
of form factors). Provided ()3.32j) holds, each term of the series has to be finite, what implies an 
upper bound on the asymptotic behaviour of the form factors. One can check that chosing the 
value of p yielding the most stringent bound gives the final result 

lim F*...^,... A^exp^l) (3.33) 

|0i|— >oo 

with 

Y* < ^ . (3.34) 

Explicit expressions for all the form factors of a given operator can be obtained in many cases. 
The problem of summing series like (j3.22|) . however, remains unsolved for interacting theories. 

6 Any operator with (Ol^l^-a) 7^ can be taken as creating operator. In order to get rid of spin factors, we 
always choose <j> a to be scalar. 

7 We denote by {■ ■ -) c the connected correlators. 
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Due to the exponential dumping factor, (|3.22|) is a large distance expansion. A complementary 
short distance expansion is provided by perturbation theory around the conformal point |26j . 
Direct comparison shows that both expansions converge very rapidly and that the first few 
terms provide a very good matching at the intermediate scales where the crossover from the 
short distance power law behaviour to the large distance exponetial decay takes place. 

Integrated correlators like (|3.31j) are what one needs to compute in many physical applica- 
tions (see e.g. section 5). In these cases the contribution of short distances is suppressed by 
the powers of \x\ in the integrand, and the form factor expansion alone is sufficient to obtain 
precise estimates. A check of the convergence can be done through sum rules allowing to recover 
conformal data from the off-critical correlation functions. Conformal theories in two dimensions 
are labelled by a number C called 'central charge' 7 j. It is a consequence of Zamolodchikov's C- 
theorem |27| that the central charge of the conformal limit of a massive theory can be expressed 



where Q(x) is the trace of the energy- momentum tensor. Similarly, the scaling dimension of a 
relevant operator <J>(a?) is computable through the sum rule [25] 



4 Integrable directions in the scaling Ising model 
4.1 Zero magnetic field 

As we have already seen, the scaling limit of the Ising model without magnetic field is described 
by the action ()2.12l) with h = 0, namely by a theory of free neutral fermions with mass m 
proportional to \t\. Depending on the sign of r, this theory describes the ordered or disordered 
phase. To be definite, we will work in the high-temperature phase (r > 0) and will use the 
duality transformations e <-> — e and a *-* [i to obtain the results in the ordered phase. 

The free fermionic theory provides a particularly simple example of integrable field theory. 
It contains a single neutral particle A with fermionic statistics and S'-matrix equal to 1. In 
order to fit the conventions on form factors of the previous section, however, we will think of the 
particle as created by a bosonic operator and absorb the anticommutativity into the 5-matrix, 
so that the two-particle scattering amplitude reads 



In the high-temperature phase we are considering, the particles correspond to local spin excita- 
tions and the spin operator a is naturally identified as the scalar creation operator 8 . 

8 In the low-temperature phase the excitations are kinks interpolating between the two degenerate vacua. Such 
topologic excitations are non-local in the spin degrees of freedom and the creation operator corresponding to the 
amplitude jUJ is the disorder operator p. 



as [23] 




(3.35) 




(3.36) 



S = -1 . 



(4.1) 
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For the present case in which we have only one type of particle the notation for form factors 
can be simplified to 

F*(0!, ...,9 n ) = (O|*(O)|A(0!) . . . A{6 n )) , (4.2) 
so that the form factor equations become 

F*{9 X + A, . . . ,9 n + A) = e sA F*(0i, ...,e n ) (4.3) 

F®(0x, . . . ,9i, 9 i+ i, . . . ,9 n ) = —F®(0i, . . . ,9 i+ i,9i, . . . , 9 n ) (4.4) 

F*{9 X + 2ivr, 9 2 ,...,9 n ) = h, a F*(9 2 , ■■■,9 n ,9 1 ) (4.5) 

Re Se >=e+i* F* +2 {9', 9,9 1 ,...,9 n )=i[l- (-l) n h, a ] F*{9 1: ...,9 n ) (4.6) 



(equation (j3,27|) plays no role in absence of bound states). We know from the discussion of 
section 2 that 

The solutions [HHI 

F®(9 1 ,...,e n ) = c 1 mF*(9 1 ,...,9 n ) = -2inm 2 5 n>2 sinh^-^ (4.8) 
F 2 Va(0i, • • • , 2n+1 ) = i n F? J] tanh (4.9) 

i<j 

F»Jf)x, • • • , M = II tanh ( 4 - 10 ) 

are the only ones satisfying the asymptotic bound ()3.34j) for the case of relevant scalar operators. 
This is a simple example of how the form factor approach reveals the operator content behind 
a scattering theory. A counting of solutions including descendents confirms the correspondence 
with the operator space of the conformal point I n writing (|4.8j) we used the fact that, 

being the operator which breaks conformal invariance, e is proportional to the trace of the 
energy-momentum tensor (ci is a dimensionless constant) for which the normalisation condition 

F®(9 + i7r,9) = 27rml (4.11) 

holds. The vanishing of the F£ with n / 2 matches the fact that the energy operator is bilinear 
in the free fermions. The results for a and \x are non-trivial and reflect the non-locality of these 
operators with respect to the fermions. The form factors of a with an even number of particles 
have to vanish because the particles are odd under a — > —a; the invariance of fi under this 
symmetry induces the vanishing of the FJ^, i- 

It follows from (|4.8j) that the form factor expansion of a two-point function involving the en- 
ergy operator contains only one term. In particular, the sum rules (|3.35J) and (|3.36|) immediately 
give the exact results 9 (Fq = (fi)) 

„ 3 f°° Ja sinh 2 1 

C = - d9 -r- = - (4.12) 

2 Jo cosh 4 # 2 v ; 



3 The integral in U.S. 36^ diverges for $ = s. See [29] on this point. 
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x„ 



1 

2^ 



00 „ sinh 2 9 
ad 



(4.13) 

cosh 3 8 v ; 

Form factor series with an infinite number of terms are instead obtained for the correlators 
(era) and (pp). Due to the particularly simple form of the matrix elements ()4.9j) and (|4.1U|I 
these series can be resummed [221 ESI- We quote here the final result of this procedure which, of 
course, reproduces that originally obtained in |2] from the lattice solution. The correlators can 
be written in the form 

( ( r(xV(0)) = ( / u) 2 e- E Wsinhi X (t/2) 



(M^)M0)) = (/i) 2 e- S Wcosh^(t/2), 



where t = m\x\ and 



m = \ 



t/2 



dpp {d P xf ~ 4 sinh 2 x 



The function x ls the solution of the differential equation 10 

1 

— i 

P 

satisfying the asymptotic conditions 



d 2 p X + -d P X = 2 sinh 2% 



x(p) 



In p + constant . 



X (p) ^ -K (2p) , 



7T 



P 



OO . 



Both (gHU) and (|4~T^|) behave as 



Ci 



Ci 



\2X a 



(4.14) 
(4.15) 

(4.16) 

(4.17) 

(4.18) 
(4.19) 

(4.20) 



as \x\ — * 0. When \x\ — > 00 the correlators decay exponentially: (era) vanish while (pp) ap- 
proaches (p) 2 - If the operators are normalised in such a way that Cj = 1, then 



(p) = ±m 1 / 8 2 1 / 12 e- 1 /^3/2 ( 

A being the Glaisher constant 

^ = 2 7 / 36 vr- 1 / 6 exp 



1 2 
3 + 3 Jo 



dx lnT(l + x) 



1.282427129. 



(4.21) 



(4.22) 



These results for the correlation functions of the scaling Ising model without magnetic field 
were also derived in [34, using the theory of monodromy preserving deformations of ordinary 
differential equations. Very recently a simpler derivation based on the Ward identities associated 
to the conservation laws of the free fermionic theory has been given in [33] . The convergence of 
the form factor series was analysed numerically in j2SJ. A study of four-point functions in the 
form factor approach can be found in |36j . 

10 Equation 14.171 becomes a Painleve III equation for the function r\ — e~ x . 
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4.2 Non-zero magnetic field at T = T c 

As anticipated in section 2, the Ising field theory 1)2.12)) with h ^ is integrable only when r = 0. 
The counting argument applied to this purely magnetic perturbation of the Ising conformal point 
shows that conserved quantities of the form 1)3. 4 Jl exist for spin s = 1,7, 11, 13, 17, 19 [""]. These 
are expected to be the first representatives of the infinite set 

s = 1,7,11,13,17, 19,23,29 (mod 30) . (4.23) 

It can be noted that this set of numbers coincides with the exponents of the Lie algebra E% 
repeated modulo the Coxeter number of the algebra. The spectrum of conserved spins ()4.23)) 
and its relation with the algebra Eg were first predicted by V. Fateev [3*71 13*5] . 

Since the Ising field theory with non-zero magnetic field does not possess any internal symme- 
try, its mass spectrum is guaranteed to be non-degenerate. Hence we know that the two-particle 
scattering amplidutes of the integrable theory must be of the form ()3.17)) . A. Zamolodchikov 
looked for the minimal solution 11 of the bootstrap equations ("3.21)1 satisfying the conservation 
equations 1)3.20)1 with the set of spin values ()4.23)) [01 . He found that this requirement leads to 
a bootstrap system closing on 8 particles A a (a = 1, ... ,8) with masses 



m 2 


7T 

= 2m\ cos — = 
5 


= (1.6180339887..) mi 


m 3 


7T 

= 2m\ cos — 


= (1.9890437907..) mi 


m± 


7vr 

= 2mo cos — 
30 


= (2.4048671724..) mi 


m 5 


2vr 

= 2mo cos — 
15 


= (2.9562952015..) mi 


m 6 


7T 

= 2mo cos — 
30 


= (3.2183404585..) mi 


m-j 


= 4m 2 cos - cos — = (3.8911568233..) mi 
5 30 V ' 

7T 27T 

= 4m 2 cos - cos — = (4.7833861168..) m x 
5 15 


m 8 



in units of the lightest mass m%. The scattering amplitudes can be written as 

s ab (0)= n (w^r (424) 

in terms of the building blocks ()3.18)k the complete list of indices 7 and exponents p 7 is given 
in Table 1. For example, the ^4i^4i and A\A 2 amplitudes read 

M0)=*2/3(0)M )*1/15(0) ( 425 ) 

Su{0) = hM0)t 3/5 (9)t 7/15 (9)t 4/15 (9) . (4.26) 

Equation 1)3. 14|) shows that the particles A\, A 2 , A3 {A\, A 2 , A3, A4) appear as bound states 
in the ^4i^4i {A\A 2 ) scattering channel. 

11 In general, the solution of the S-matrix equations of section 3.2 is not unique. There exists, however, a 
minimal solution possessing the smallest number of zeros and poles in the physical strip Imt9 £ (0, ir). 
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3 

(2) 

3 4 

(14) (8) 

(13) (3) (H) 2 
(17) (A) (7) (15) 

(14) (4) (10) 2 (12) 2 
(9) (7) 2 (13) 2 (15) 

(5) (9) 2 (11) 2 (13) 2 (15) 
3 (6) 2 (8) 2 (10) 2 (12) 2 
(14) (8) (2) (12) 2 
(9) (7) 2 (13) 2 (15) 
(5) (9) 2 (11) 2 (13) 2 (15) 
3 (6) 2 (8) 2 (10) 2 (12) 2 
(19) 3 (13) 3 (3) (7) 2 (9) 2 (15) 
3 (17) 3 (11) 3 (5) 2 (7) 2 (15) 2 
3 (4) 2 (6) 2 (10) 4 (12) 4 (16) 4 
3 (14) (12) 3 (4) (2) 2 
3 (6) 2 (8) 2 (10) 2 (12) 2 
(2^) 3 (13) 3 (5) (3) 2 (11) 4 (15) 



(18) 3 (8 8 ) 3 (10) 2 (16) 4 
(4) 2 (6) 2 (10) 4 (12) 4 (16) 4 



(17) 5 (7) 4 (9) 4 (11) 2 (15) 3 
(16) 3 (12) 3 (2) (6) 2 (8) 2 
(19) 3 (9) 3 (5) 2 (13) 4 (15) 2 
(4) 2 (6) 2 (10) 4 (12) 4 (16) 4 

(18) 5 (14) 5 (4) 2 (8) 4 (10) 4 
(21) 5 (3) 2 (7) 4 (11) 6 (13) 6 (15) J 
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8 8 



(22) 3 (20) 5 (12) 5 (2) 2 (4) 2 (6) 2 (16) 4 
(27) (25) 3 (17) 5 (7) 4 (9) 4 (11) 4 (15) 3 
(29) (25) 3 (21) 5 (3) 2 (7) 4 (11) 6 (13) 6 (15) 3 

3 4 5 8 

\3 fiA\5 M S^7 / a \6/ 



(28) (26) 3 (24) 5 (18) 7 (8) b (10) b (16) 8 
(24) 3 (20) 5 (14) 5 (2) 2 (4) 2 (8) 4 (12) 6 

(28) (26) 3 (22) 5 (16) 7 (6) 4 (10) 6 (12) 6 

(29) (27) 3 (23) 5 (21) 7 (5) 4 (11) 8 (13) 8 (15) 4 
(26) 3 (24) 5 (20) 7 (2) 2 (8) 6 (12) 8 (16) 8 

(29) (27) 3 (25) 5 (23) 7 (19) 9 (9) 8 (13) 10 (15) 5 
(28) 3 (26) 5 (24) 7 (22) 9 (20) 11 (12) 12 (16) 12 



Pi 



Table 1: 5-matrix of the scaling Ising model in a magnetic field at T = T c . A factor (i 7 /3o(#) 
in the amplitude S ab {9) corresponds to each term (7) P7 (p 7 = 1 is omitted; the building blocks 
t a {6) are defined in (13,181) ). The superscript c above (7) indicates that the pole at 9 = iir^/SO 
in the amplitude S ab (9) corresponds to the particle A c appearing as bound state in the A a A b 
scattering channel. 



Most of the amplitudes in Table 1 contain higher order poles that we did not discuss in the 
previous section. Within the framework of the analytic 5-matrix, each singularity is expected 
to have a physical interpretation. The higher order poles of the scattering amplitudes were 
recognised in jHS] (see also [313 |^ ) as the singularities that in two dimensions are associated 
to processes in which more than one particle propagates in the intermediate state 12 . Figures 7 
and 8 show the processes of this kind leading to second and third order poles at 9 = i(p, 

V = <d + u e db -TT, (4.27) 
in the amplitude S ab (9). The angle 77 in Figure 7a is 

r7 = 7r-t&-t&€[0,7r), (4-28) 

and irj is the rapidity difference between the intermediate propagating particles A c and A e ; the 
diagram of Figure 7b corresponds to the limiting case r/ = 0. One has 

S ab (9 ~ up) ~ ^^ffi^ (4-29) 

pa pb pfo pa 

S ab (9 ~ Up) ~ y_ e ,^ )2 /e (4.30) 



2 In four dimensions these processes give rise to anomalous thresholds rather than poles |21) . 
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Figure 7: Scattering patterns associated to the second order poles of the S-matrix. 



for the diagrams of Figures 7a and 7b, respectively, and 

(ya yb yf \2 

Soh (9 * up) * i { (4.31) 

for the direct channel third order poles of Figure 8. More generally, a pole of order P — 2L 
corresponds to a diagram with P internal lines and L loops |4()j . 

We said that the S-matrix of Table 1 is the minimal solution to the constraints discussed so 
far. Although it is extremely natural to expect that this minimal S'-matrix is the one describing 
the magnetic Ising model at T = T c , the conjecture needs to be checked. Al. Zamolodchikov 
showed that the conformal limit of an integrable field theory can be identified using the S'-matrix 
as the only imput. The method, known as thermodynamic Bethe ansatz, allows the computation 
of the central charge of the conformal limit through the study of the thermodynamics on a 
cylindrical geometry When applied to the S-matrix of Table 1 it yields the expected result 
C = 1/2 jug. 
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Figure 8: Third order pole of the scattering amplitudes. 



A different way of confirming that the S-matrix of Table 1 does correspond to the Ising model 
is that of showing that it describes an integrable quantum field theory containing two relevant 
scalar operators besides the identity 13 . Addressing this problem means taking Zamolodchikov's 
S'-matrix as the imput of the form factor equations (|3.25M3~2U|) and looking for scalar solutions 
which behave asymptotically as (j3.33j) with Y$ < 1 [H1E3] . The non-degenerate spectrum and 
the "(/3 3 -property" ^ show that the theory possesses no internal symmetries, so that all 
the form factors of scalar primaries other than the identity are expected to be non- vanishing. 
For operators local with respect to the particles (l^^ a = 1), the general solution of equations 
(|M.2fi|) and (|3.28[) with the required pole structure can be written in the form 



rpmin I n 

<i (coah(^))^ D aiaj 



*S_«.(*i, ...,9 n ) = Qt..^, . . . A) II - 3 ■ (4-32) 



Here F™ n {9) is a solution of the equations 

F{9) = S ab (9)F(-9) (4.33) 

F(9 + 2m) = F{-9) (4.34) 

free of zeroes and poles for Im.9 G (0, 2tt). In the denominator, the factors cosh ( ' 2 j J introduce 
the annihilation poles, while D a b(9) takes care of the dynamical poles in the A a Ab scattering 
channel through factors of the type (cosh# — cosu^ b ) (see later for the precise form). Finally, 
the QaL^on are entire functions of the rapidities, invariant under exchanges 9\ <-> 9j and (up to a 
13 This requirement uniquely identifies the Ising model among the theories satisfying reflection positivity. 
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sign) 27ri-periodic in all Ofs. They are subject to (|3.25f) and to the residue equations Q3.27JI and 
([3.29ft which relate functions with different n. These equations, however, hold for any operator 
with a given spin, so that further constraints are needed to select specific solutions. We will see 
in a moment the role played in this respect by the asymptotic bound ()3.34|) . 

The form factors are determined starting with the first non-trivial case (the two-particle 
one) and then using the residue equations to fix the matrix elements with an higher number 
of particles (form factor bootstrap). The specialisation of ()4.32[) to n = 2 and scalar operators 
reads 

F* b (8) = jf^F% in (0), (4-35) 

where we made the identifications F%{0i,0 2 ) = F^(9i - 9 2 ), QfU^A) = Q%,{0\ - #2), and 
took into account the vanishing residue on the annihilation pole in the two-particle case when 
h,4> a = !• 

The functions F^ m (9) with the properties specified above can be written as 



ab 



\6)= -*sinh- J] (jVaoWr ' ( 436 ) 



where 



oo dt cosh(a-ijt „(m-6)t. 



solves the equations 



T a {9) = -t a (8)T a (-0) (4.38) 
T a {d + 2m) = T a {-e) ; (4.39) 



the property 



S ab (0) = {-lf ab (4.40) 

of the S-matrix originates the factor in front of the product in (j4.3fij) . 

As we said, the denominator of (|4.35j) takes care of the dynamical poles. While equation 
(|3.27|) says how to deal with the simple poles of the S-matrix, the case of higher order poles 
needs to be discussed Due to crossing symmetry, poles appear in pairs in the scattering 
amplitudes. For poles of even order both residues are positive and there is no way to distinguish 
between a direct and a crossed channel. Then a diagram of the type shown in Figure 7 must 
exist for each double pole at 9 = iip in the amplitude S a b(9). In the vicinity of such a pole, the 
form factor behaves as 14 (see Figure 9) 

F%(9 ci<p)ci ^ASce^F?^) = . rjAFtiiv) 
a ' ' 9 — iip 9 — iip 

(this result also holds for 77 = 0). 

14 The pole in the form factor is simple because the diagram of Figure 9 has a single triangular loop. 
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Figure 9: Form factor singularity associated to a double pole of the S-matrix. 




Figure 10: Form factor singularity associated to a direct channel triple pole of the 5-matrix. 



A third order pole in the scattering amplitude can be seen as originating from a double pole 
when r\ = u[ e . Then the direct channel singularity at 9 = i<p in the form factor is obtained using 
(|3~27|) into (|OT|) and reads (Figure 10) 

F%(6 * up) * - T ^f}f^ F f > ( 4 - 42 ) 

while the crossed channel pole at 9 = i(ir — <p) remains simple. 

The analysis of poles of higher order can be done along similar lines and leads to the expres- 
sion 

Dab(0) = I] (?W#)) t7 (Vi-WsoW)' 7 , (4-43) 

where 



cos not — cosh 9 

2 



and 

z 7 = n + 1 , j 1 = n , if p 7 = 2n + 1 



(4-45) 

ij = n , j-y = n , if p 7 = In . 
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We are now in the position of dealing with the functions Q%,(0), the only piece of (|4.35j) 
which carries information about the operators. As functions of the rapidity difference, they 
must be even, 2-7ri-periodic, free of poles and exponentially bounded. Hence we can write them 
as 

AT* 

Qtb(8) = E c % coshfe 6 • (4-46) 

k=0 

The condition 

= (4-47) 

follows from Q3.26JI and S* b (6) = S a b(—6), and ensures that the coefficients cj^$ are real. These 
coefficients are the only unknowns we are left with. 

Let us look for solutions corresponding to relevant operators. This means that the asymptotic 
behaviour (|3.33|) has to hold with Y$ < 1. Since 

T a {0) ~ exp(|0|/2) , |0|->oo, (4.48) 

it is straightforward to check on (|4.35|) that, in particular, 

Nf x < 1 . (4.49) 

Hence, the initial condition of the form factor bootstrap for relevant scalar operators allows for 
two free parameters, i.e. the coefficients cfj^ and c^'^,. It can be checked that the number of 
free parameters does not increase when implementing the bootstrap. For example, since 

N? 2 <2, (4.50) 

considering Ff 2 {6) brings in three more coefficients. On the other hand, the amplitudes S\\{9) 
and S\2{0) possess three common bound states, what yields the three equations 

-^Res^^F*^) = -LRes e=i ^ 2 if 2 (#) , c = 1, 2, 3 (4.51) 

1 11 1 12 



(k) (k) 

which determine the three c 12 $ in terms of the two cii Going on and using also the conditions 
on higher order poles, a number of residue equations larger than the number of new coefficients 
is available in many cases. It turns out, however, that the extra constraints are automatically 
fulfilled so that the two initial parameters propagate untouched through the bootstrap procedure 

puni. 

Since all the form factor equations used so far are linear in the operator <I>, the interpretation 
of this result is simple: the Zamolodchikov's S"-matrix describes an integrable field theory with 
two scaling relevant scalar operators and $2 (plus the identity) and, up to an additive 
constant, the operators $ selected up to now can be written as 

$ = a$i + /3$ 2 (4.52) 

with a and (3 real parameters. We already said that this condition uniquely identifies the Ising 
field theory (|2.12|) in which the two scaling operators are a and e; the known results about the 
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integrable directions imply r = 0, h ^ 0. Hence, the analysis of the operator content confirms 
that the Zamolodchikov's S"-matrix corresponds to the integrable direction of the magnetic Ising 
model. 

Obviously, the next issue is that of identifying the solutions corresponding to the scaling 
operators a and e. Having already exploited all the constraints on form factors discussed so 
far, it is clear that this task requires some new physical input. Since the scaling operators 
are characterised by their behaviour close to criticality, it is natural to look once again at the 
asymptotic properties of form factors. More precisely, we consider the limit j^H] 

JSR. F* Us ,..„ (ft + 2, . . . , ft + f ,ft - ° . . . , «i - £) , (4.53) 

where 

*fc=/|\> k = M (4.54) 

are the scaling operators entering (|4.52j) rescaled by their vacuum expectation value. It follows 
from (j3.10j) that shifting a rapidity 9 by ±a/2 and rescaling the mass as m a = M a e~ a / 2 produces, 
when a — ► +oo, a massless particle with energy p° = ip 1 = H^e^ e (right- or left-mover, 
depending on the sign of momentum). Hence, (|4.53|) corresponds to the massless limit towards 
the conformal point in which the first r particles become right-movers and the remaining I 
particles left-movers. The property 

lim S ab {9 + a) = 1 (4.55) 

of the scattering amplitudes of Table 1 shows that a right-mover and a left-mover do not interact 
in such a conformal limit. We conclude that the limit 1)4. 53|) produces a factorisation into two 
massless form factors, one with r right-movers and one with I left-movers. On the other hand, 
a massless form factor with all right (left) movers is obtained from (|4.53|) with 1 = (r = 0); in 
such a case all rapidities are shifted by the same amount, and (|3.25|) shows that this massless 
form factor actually coincides with the massive one. Then we conclude that the asymptotic 
factorisation property 

lim l-ti. (9x + a,...,e r + a,9' 1 ,...,9' l ) = ^ fc ..a r • • • , k • • • , e[) (4.56) 

(X — >"-pOO x t 

holds in the original massive theory. Notice that the particular cases r = and/or I = 
require the normalisation Q4.54JI . The factorisation argument leading to (|4.56|) applies to scaling 
operators and not to linear combinations like (|4.52|) mixing operators with different scaling 
dimensions: in the latter case one of the coefficient a, (3 is dimensionful and vanishes in the 
massless limit. 

We are now in the position of identifying the initial conditions of the form factor bootstrap 
corresponding to the scaling operators $i and <I>2- This amounts to finding the values of the 
ratio 

(0) 

z* = (4.57) 
c ii,$ 
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for which the factorisation conditions (|4.56j) are fulfilled (this quantity is universal because does 
not depend on the normalisation of the operator). The equation 

4r I™ F^(9) = -L. lim F* k (9) (4.58) 

is a consequence of (|4.56j) and gives a quadratic equation for z§ whose solutions are pOJ 

= 4.869840.. 
z$ 2 = 1.255585.. . 

Being the operator which perturbs the conformal point, a is proportional to the trace of the 
energy-momentum tensor. It is not difficult to see that the conservation of the latter induces 
the factorisation of a kinematical term in the form factors of a. In particular, the two-particle 
form factors F° h (6) contain the factor [5] 

(2 , 2\ 1_ *a!> 

cosh0 + ^±^) . (4.59) 

2m a m b J 

It can be checked that the form factors originating from the initial condition have this 
property. Then we conclude 

$i = a 
$2 = e . 

Having determined the initial conditions of the form factor bootstrap for the two relevant 
operators, all their form factors can in principle be computed using the residue equations on 
dynamical and kinematical poles. The results for several two-particle form factors are given in 
Table 2; Table 3 contains the full list of one-particle matrix elements Referring to the 

normalisation- independent ratios (|4.54jl . all these numbers are universal. One can check that 
the factorisation conditions prescribed by (|4.56|) are satisfied. Of course, the bootstrap can be 
continued to determine the remaining two-particle form factors as well as those with more than 
two particles (for example, F[ u (9i, 62, 63) is computed in [S]). 



The content of Tables 2 and 3 is sufficient to develop the large distance expansion (|3.22|) of 
two-point correlators including all terms of order lower than e - 3m iM. J n particular, the first few 
terms are 

($ fe (x)d,(0)) = 1 + - E F^FpK {m a \x\) + 0(e~ 2m ^) , (4.60) 

where Kq(z) is a Bessel function. 

The sum rules (|3.35[) and ()3.36|) can be used to test the convergence of the form factor 
expansion. We recall that O ~ a and that the normalisation of the form factors of O is fixed by 
the condition (|4.11|) . The asymptotics 

(a(x)a(O)) ~ ^ 
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a 




11 

c ll 


-2.093102832 
-10.19307727 


-70.00917205 
-87.90247670 


r 2 
c 12 

r 1 
r° 

1Z 


-7.979022182 
-71.79206351 
-70.29218939 


-466.3008246 
-1307.331521 
-853.2803886 


c 13 

r 2 
c 13 

c 13 

lo 


-582.2557366 
-6944.416956 
-13406.48877 
-7049.622303 


-43021.45153 
-182413.2733 
-241929.7678 
-102574.1349 


r 2 
c 22 

c 22 
r° 

ZZ 


-21.48559881 
-333.8125724 
-791.3745549 
-500.2535896 


-2193.896354 
-10870.05277 
-16161.44508 
-7510.235388 


r 3 
c 14 

r 2 

c 14 
c 14 


22.57778351 
318.7122159 
672.2210098 
377.4586311 


2074.636471 
9881.413381 
14357.04570 
6568.762583 


cf, 

15 

c 15 
r 2 

c 15 


-260.7643072 
-4719.877128 
-15172.07643 
-17428.22924 
-6716.787925 


-30333.56619 
-198757.2340 
-447504.5720 
-422808.9295 
-143743.2050 


( A 
23 

r 3 
c 23 

r 2 
c 23 

c 23 

c 23 


-92.73452350 
-1846.579035 
-6618.297073 
-8436.850082 
-3579.556465 


-11971.94909 
-81253.72269 
-186593.8661 
-178494.3378 
-61194.62416 


J> 
c 33 

r 4 

LOO 

oo 

c 33 

r 2 
c 33 

c 33 

c 33 


-1197.056497 
-30166.99117 
-150512.4122 
-301093.9432 
-267341.1276 
-87821.70785 


-195385.7662 
-1743171.802 
-5603957.324 
-8422606.859 
-6035102.896 
-1668721.004 
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c 25 


1425.995027 


289831.4882 


c 25 


44219.03877 


3275586.983 


c 25 


286184.1535 


13872077.63 


c 25 


788413.2178 


29236961.96 


c 25 


1078996.488 


32979257.31 


c 25 


725356.4417 


19100224.04 


c 25 


191383.5734 


4471623.121 


c 17 


190.8548023 


30394.23374 


r 4 
c 17 


4633.706068 


274294.8033 


C?7 


21406.72691 


897781.3229 


C?7 


39514.82959 


1375919.456 


c 17 


32456.91939 


1004969.466 


c 17 


9906.265607 


282938.1974 


C 44 


-7249.785565 


-1830120.693 




-276406.7236 


-25699492.93 


C44 


-2299573.212 


-138411873.8 


C 44 


-849276.3526 


-384776478.8 


r 3 
C 44 


-16615618.39 


-608371427.1 


C44 


-17950817.11 


-553818699.0 


C44 


-10139089.36 


-270964337.7 


C 44 


-2341590.241 


-55283137.91 



Table 2: Coefficients of some of the polynomials (|4.46jl for the primary operators a and i 
(parenthesis on the upper index of the coefficients are omitted; = <3?/ (<£)). 





a 


£ 


Ft 


-0.64090211.. 


-3.70658437.. 


F 2 


0.33867436.. 


3.42228876.. 


F3 


-0.18662854.. 


-2.38433446.. 


F* 


0.14277176.. 


2.26840624.. 


F 5 


0.06032607.. 


1.21338371.. 


F 6 


-0.04338937.. 


-0.96176431.. 


F 7 


0.01642569.. 


0.45230320.. 


F 8 


-0.00303607.. 


-0.10584899.. 



Table 3: The one-particle form factors of the primary operators a and i (3> = 
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Ci 


0.472038282 


c 2 


0.019231268 


c 3 


0.002557246 




0.003919717 




0.000700348 


Cl2 


0.000974265 


c 5 


0.000054754 


Cl3 


0.000154186 


Cpartial 


0.499630066 



Table 4: Central charge from a partial sum of the form factor expansion of the correlation 
function in the sum rule ()3.35j) . Cab,, denotes the contribution of the state A a A .. . The exact 
result is C = 1/2. 

Mx)e(0)) ~ t^T (ff) (4-61) 
| a; | 

(8(x)e(0)) ~ ^§ 
\x\- 

as | a; | — > follow from (|2.4|) and (|2.7j) and are helpful in the interpretation of the results of 
Table 4 and 5 for the sum rules. The faster convergence of the central charge sum rule is 
expected since the integration in |x| 3 d|a;| strongly suppresses the contribution of short distances. 
In the scaling dimension sum rule the suppression is lower by a factor \x\ 2 and more effective 
for X a than for X £ since (era) is less singular than (ere). 

It is of obvious interest to compare the results of integrable field theory discussed in this 
section with the numerical results for the original lattice model (|2.1jl with T = T c . Actually, 
the first numerical investigation was performed in |U] on the Ising quantum spin chain with the 
purpose of confirming the Zamolodchikov's mass spectrum, what was done with good precision 
for the first few masses 15 . Concerning correlation functions, after the early studies of (481 149j . the 
most accurate Monte Carlo simulations have been performed in jSJJ. In this work the numerical 
data have been used to test the form factor expansion as well as the corrections to the short 
distance behaviour computed in the analysis has confirmed the good convergency properties 
of both expansions. 

Caselle and Hasenbusch used a numerical diagonalisation of the transfer matrix of the lattice 

model to evaluate the first few one-particle form factors 52: . The comparison of their results 

(Table 6) with those of Table 3 provides a direct check of the form factor computations reviewed 

in this section. The transfer matrix method has also provided a test of two-particle form factors. 

15 The full spectrum has been recovered in |45| from the exact solution of the RSOS model of 0, which is in 
the same universality class of the magnetic Ising model. The spectrum of the field theory 12.121 with r = was 
also studied in |46| by numerical diagonalisation of the Hamiltonian on the (suitably truncated) space of states 
of the conformal theory |47) . 
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a 


£ 


Ai 


0.0507107 


0.2932796 


A 2 


0.0054088 


0.0546562 


A 3 


0.0010868 


0.0138858 


A n 


0.0025274 


0.0425125 


A 4 


0.0004351 


0.0069134 


A12 


0.0010446 


0.0245129 


A 5 


0.0000514 


0.0010340 


A13 


0.0002283 


0.0065067 


Apartial 


0.0614934 


0.4433015 



Table 5: Conformal dimensions A$ = X$/2 from a partial sum of the form factor expansion 
of the correlation functions in the sum rule (|3.36j) . A a &.. denotes the contribution of the state 
A a A b .. . The exact results are A CT = 1/16 = 0.0625 and A £ = 1/2. 



I*i| 

1*21 



a 



0.6408(3) 
0.325(25) 



3.707(7) 
3.38(7) 



Table 6: Lattice estimates of one-particle form factors obtained in |52j . 



Indeed, one-point functions on a cylinder of circumference R behave as 
where the amplitudes 



A $ ^ F aa( i7r ) 



R=oo 



are given in Table 7. The avalaible numerical estimates are |53| 

-8.11(2), A\ = -17.5(5) 



R — > oo 



(4.62) 



(4.63) 



(4.64) 





a 


e 


A x 


-8.0999744.. 


-17.893304.. 


A 2 


-21.206008.. 


-24.946727.. 


A 3 


-32.045891.. 


-53.679951.. 



Table 7: Universal amplitudes ruling the leading finite size corrections of one-point functions 
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5 Ising universality class 



Two statistical mechanical systems may differ for a number of microscopic features (e.g the 
type of lattice or the number of neighbouring sites entering the interaction term) and still be 
characterised by the same internal symmetry. In this case, the systems will exhibit the same 
critical behaviour nearby a second order phase transition point associated to the spontaneous 
breaking of the symmetry. As it is said, they belong to the same universality class. Quantum 
field theory deals directly with the continuum limit in which the microscopic details become 
immaterial, and is the natural framework for describing universality classes. In particular, 
the action (|2.12[) describes the universality class of the two-dimensional Ising model, and all 
the results discussed so far with reference to this action are universal. Traditionally, however, 
the quantitative characterisation of universality in statistical mechanics is made studying the 
behaviour of the thermodynamical observables in the vicinity of the critical point. In this section 
we show how integrable field theory allows to complete the list of canonical universal quantities 
of the Ising universality class • 

The usual procedure (see |56j and references therein) is that of introducing critical exponents 
and critical amplitudes through the relations 

'f±\r\- u , t^0±, h = 

< 

. fc\h\- v ; r = 0, h^O 

( (A±/a) |r|" a , r^0 ± , h = 
[ (A c /a c )\h\-^, r = 0, h^O 

( B{-rf, t-»0" h = 
((\h\/D)^ s , r = 0, h^O 



£ = 



C = 



\M\ = 



fr±|r|-T, r^0 ± , h = 

X = 

{ r c |/i|->, r = 0, h^O 

where £ is the correlation length, C the specific heat, M the magnetisation and x the suscep- 
tibility; r and h are the deviation from critical temperature and the magnetic field entering 



the action (|2.12|) . The observables above are related to the one- and two-point functions of the 
scaling operators a and e. After introducing the free energy per unit area 

f = ~lnZ, (5.1) 
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we have 16 

C = "0= / d 2 x(e(x)e(0)) c (5.2) 

M = -% = (*) (5.3) 

x = ~% = /^K*MO))c- (5.4) 

As for the correlation length, it is common to distinguish between the 'true' correlation length 
£t defined as 

lim {a(x)a(0)) c ~ e" 1 ^ , (5.5) 

— s-oo 

and the second moment correlation length 

&nd = (J- J d 2 x\x\ 2 (a(x)a(0)) c y 2 . (5.6) 

We will distinguish the critical amplitudes for the two types of correlation length through the 
superscripts t and 2nd. 

If M ~ l/£ is a mass scale, we know that / ~ M 2 , r ~ M 2 ~ x ^ and h ~ M 2_Xct . Then, 
comparison of (|5.2jl - (|5.f)jl with the definitions of the critical exponents gives 

v = 1/(2 -X e ) = l 

v c = 1/(2 - X a ) = 8/15 

a = 2(1 -X e )v = 

a c = 2(1 - X £ ) v c = 

P = X a v = l/8 

5 = l/(X a v c ) = 15 

7 = 2(1 - X CT ) = 7/4 

7c = 2(1 - X CT ) v c = 14/15 . 

These well known results express the fact that the critical exponents are universal and determined 
by the scaling dimensions of the operators a and e. They also account for the usual scaling and 
hyperscaling relations 

a + 2/3 + 7 = 2 (5.7) 

a + 2v = 2 (5.8) 

7 = P(5-1) (5.9) 

a c = a / (38 (5.10) 

v e = u/fiS (5.11) 

7c = 1 " 1/5 • (5-12) 



16 We recall that {■ ■ -) c denotes connected correlators. 



32 



Contrary to the exponents, the critical amplitudes are not determined by conformal field 
theory and are not themselves universal. Roughly speaking, they depend on the scales we use 
to measure the temperature and the magnetic field (metric factors). It is possible, however, 
to combine different amplitudes in such a way that the metric factors cancel leaving universal 
quantities. Such universal amplitude combinations provide the canonical characterisation of the 
universality of the scaling region surrounding the critical point. 

A simple way of canceling the metric factors is to consider ratios like r + /r_, fc nd /fc an d 
similar. More sophisticated cancelations are obtained exploiting the relations (|5.7[) - (|5.12|) among 
the critical exponents. So, to (|5.7|) - (|5.11|l one associates the universal quantities [HE] 

R c = A+T+/B 2 (5.13) 

i? ? + = 4 /2 /t (5-14) 

R x = T + DB 5 - 1 (5.15) 

R A = A c D-^ +a ^B~ 2 ^ (5.16) 

Q 2 = {T + /T e ){f t JJ%V lv ■ (5-17) 
Equation ()5.12[) is instead associated to the identity 

ST C D^ S = 1 (5.18) 
which can be recovered from 1)3.36(1 with $ = a and 

8 = -2irh(2 - X a ) a , r = 0. (5.19) 

Notice that the critical amplitudes of the Ising universality class are defined along the inte- 
grable directions of the field theory (|2.12|) . As we are going to see, all of them can be computed 
exactly. We will work within the so-called 'conformal normalisation' of the operators a and e, 
which amounts to take 

CL = C[ e = 1 (5.20) 

in (|4.61j) . We know that the mass m of the particle at h = and the mass m\ of the lightest 
particle at r = can be written as 

m = C T \r\ (5.21) 
m 1= C h \h\ 8 ^ . (5.22) 

The two constants are known exactly from a modified version of the thermodynamic Bethe 
ansatz, and read EH] 

C T = 2vr (5.23) 
C h = 4.40490857.. . (5.24) 

The spontaneous magnetisation at h = was given in (|4.21|) and corresponds to the ampli- 
tude 

B = 1.70852190.. . (5.25) 
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In the magnetic direction, we can use (|4.56f) and Q5.19[) to write 

w = '- 2 ^ 2 -^'-' n^£rW T=0 (5 ' 26> 

where the form factors of coincide with those determined in section 4 for a, up to the nor- 
malisation fixed by (|4.11j) . The result for the magnetisation obtained in this way concides with 
that given by the thermodynamic Bethe ansatz |58| . and provides the amplitude 

D = 0.0253610264.. . (5.27) 

The high- and low-temperature susceptibilities at h = are obtained integrating the corre- 
lators ()4.14j) and ()4.15l) with the normalisation (|4.21|) . what gives 

T+ = 0.148001214.. (5.28) 
T_ = 0.00392642280.. . (5.29) 

The susceptibility amplitude at r = follows from (|5.18|) and (|5.27l) and reads 

T c = 0.0851721517.. . (5.30) 

It appears from (|5.5|) that the 'true' correlation lenght is just the inverse of the total mass of 
the lightest particle state entering the form factor expansion of the spin-spin correlator. Then 
we know from the analysis of section 4 that 

!l/m, t > 

h = (5.31) 
l/2m, r<0 

^ = 1/mi , r = . (5.32) 

Concerning the second moment correlation length, at h = the integration of the exact corre- 
lators gives 

r (0.99959808.. )m , r>0 
6nd = | , h = 0. (5.33) 

[ (0.31607894.. )m, r<0 
At r = we do not know the exact spin-spin correlator, but the integral entering (|5.6|) is related 
to the central charge sum rule (|3.35|) by ()5.19j) . Putting all together one finds the amplitude 



/ 2nd = \ — D 1 / 30 = 0.21045990.. . (5.34) 
V 45-7T 

Dealing finally with the specific heat, we have to take into account that, since a = a c = 0, 
the critical behaviour gets modified into 

-A ± ln|r|, h = 

(5.35) 

-A c ln\h\, r = 0. 
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= 1 

T + /r_ = 37.6936520.. 

fX/f- = 2 
fl nd /fl = 0.99959808.. 

f 2nd J j 2nd _ 3 .1 62 49504.. 

R c = 0.318569391.. 
i?+ = l/y/2n = 0.39894228.. 
find/ ft = 0.9270566.. 

R x = 6.77828502.. 
R A = 0.0250658794.. 
Q 2 = 3.23513834.. 

Table 8: The universal amplitude combinations of the Ising universality class. 



On the other hand, (|4.61j) and (|5.2()j) allow to isolate the logarithmic singularity as 

C = [ d 2 x <e(x)e(0)) c ~ 2vr f — 2vrln(Mr ) , (5.36) 

J Jro T 

where r = \x\, ro is a short distance cut-off and M is a mass scale needed to make dimensionless 
the argument of the logarithm. It will be propotional to m at ft = and to m\ at r = 0. 
Recalling Q5.21JI and (|5.22|) we conclude 

A± = 2vr (5.37) 

167T 

A c = — . (5.38) 

The values of the universal amplitude combinations which follow from these results are 
collected in Table 8. The first seven quantities involve only amplitudes computed at h = and 
are all known since the work of JSj; the last four numbers involve also amplitudes computed 
along the magnetic direction and where given in 55;. For the latter, an early lattice estimate 
coming from series expansions is R x ~ 6.78 |59j : more recently, transfer matrix techniques have 
given the results R x = 6.7782(8) and Q 2 = 3.233(4) [H21 



6 Beyond integrability 

We saw in section 3 how two-dimensional quantum field theories allow for integrable directions 
in coupling space along which a great deal of exact information can be obtained. A number 
of consequences of this fact for the scaling Ising model have been discussed in the last two 
sections. It is quite natural to ask how the exact results for the integrable directions can be 
exploited for the analytic study of the non-integrable ones [HO]. In general, we can decompose 
the action A ni of the non-integrable theory into the action Ai describing an integrable direction 
plus the contributions of the scaling operators responsible for the deviation from integrability. 
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Considering, for the sake of simplicity, the subspace spanned by only one of such operators, we 
write 

(6.1) 



Ani = Ai - A J d 2 x *(x) 



In principle, this theory can be studied by perturbation theory in A: all the corrections can 
be formally written in terms of the matrix elements of the perturbing operator \E' between 
particle states in the unperturbed (integrable) theory; as we know, these matrix elements can 
be computed exactly. 

The first order corrections are particularly simple and normally sufficient to explore small 
deviations from integrability. For example, the leading corrections to the energy spectrum and 
the scattering amplitudes for the case of neutral particles read [BIT] 



5ml = ■ 
5S ab (8) 



<*>Ia=o A + 0(A 2 



_ ■ F abab(°) 



A=0 



A + 0(A 2 ) 



A 



A=0 



sinh 6 



+ 0(X 2 ), 



(6.2) 
(6.3) 

(6.4) 



m a m b 

where £ vac is the energy density of the vacuum state and 17 

F*(in) = (A a (0)\*(O)\A a (e)) c = F*(e + in,e) (6.5) 

^L(fi - fc) = (A a (ei)A b (e 2 )mo)\M0i)A b (9 2 )) c = F* ab (e 1 + iir,e 2 + m,e 1 ,e 2 ) . (6.6) 

Going to the Ising field theory (|2.12jl . one possibility is to identify the magnetic direction 18 
7] = as the integrable unperturbed theory, and to take A^ = re as the perturbation. Then, in 
particular, the ratios 



r)=0 



(6.7) 



are completely universal and can be checked in any numerical approach giving access to the 
energy spectrum 19 . Their values for a = 1,2,3,4 follows from the results of Table 2, which give 



*22(wr) 



il=0 
»?=0 
»?=0 
»?=0 



-17.8933.. 
-24.9467.. 
-53.6799.. 
-49.3169.. . 



(6.8) 
(6.9) 
(6.10) 
(6.11) 



The first few ratios 1)6 .7|) have been checked with good accuracy by numerical diagonalisation of 

the Hamiltonian of the Ising field theory on a conformal basis of states [BU] , and in the lattice 

17 We denote by (•••[•■•[•■ -) c the connected part of a matrix element; ^ is a scalar operator. 

18 The parameter r\ which labels the renormalisation group trajectories in the plane of Figure 1 was defined in 

19 When the system is compactified on a cylinder of circumference R, the ground state energy behaves as S vac R 
for fl/C > !■ 
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model by numerical diagonalisation of the transfer matrix |61| . The absolute variations of the 
energy levels can be obtained using the result 



which refers to the conformal normalisation (|5.20[) of the operators and has been obtained in 
jH2] through quite sophisticated techniques; a previous lattice estimate agrees well with this 
exact value. 

Notice that 5m a will present an imaginary contribution at order r 2 for a = 4, ... ,8. Indeed, 
we know that at r = the masses m a with a > 3 lie above the lowest two-particle threshold 2mi. 
There is nothing wrong with this as long as the theory is integrable and all inelastic processes 
are forbidden. Moving away from r = 0, however, integrability is lost and nothing can prevent 
the decay A a — > A\A\ when a > 3. The imaginary part in the mass is inversely proportional to 
the lifetime of the unstable particle. For A4 it is given by [HOI 



with 6* determined by energy conservation; heavier particles receive similar contributions by 
other decay channels. These decay widths, as well as the corrections to the scattering amplitudes, 
can in principle be computed turning the crank of the form factor bootstrap. 

In a similar way, the region close to the thermal axis in Figure 1 can be analysed perturbing 
the h = theory with the operator \I/ = a. Let us start from the low-temperature phase 
(77 = —00). We know from l)6.3j) and Q4.1UJI that in this case the first order correction to the 
mass of the particle of the unperturbed theory is proportional to tanh(i7r/2), and then is infinite. 
This result is easily understood once we remember that the particle we are dealing with is the 
kink interpolating between the two vacua of the spontaneously broken phase at h = 0. As 
soon as we switch on the magnetic field, the degeneracy of the vacua is lifted, so that all states 
carrying a topologic charge acquire an infinite energy and are removed from the theory; the 
single particle excitations of the perturbed theory are kink-antikink bound states A n . Hence, 
the divergence of the mass correction is the signature of the confinement of the topologic charge 
and manifests itself whenever the perturbing operator ^> is mutually non-local with respect to 
the particles of the unperturbed theory 20 |6()| I64|: indeed, (|3.29j) shows that this is the condition 
for having a pole at 6 = in in F^ a (9). 

The magnetic field induces an energy density difference 



between the two vacua, v being the spontaneous magnetisation (|4.21|) of the theory at rj = —00 
(v takes the sign of h). Then, neglecting the relativistic effects, a kink and an antikink separated 
by a distance R are bounded by the potential 2vhR. The associated Schrodinger problem gives 
20 We recall that in the present case the kink is created by the disorder operator fj,. 




(6.12) 




(6.13) 




(6.14) 
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for the bound states A n the masses 



~ 2m + 



(2vh) 2 / 3 



(6.15) 



where m is the mass of the kink and z n , n = 1, 2, . . ., are positive numbers determined by the 
zeroes of the Airy function, Ai(— z n ) = 0. This non-relativistic approximation holds as long 
as m n — 2m <C m, and then gives the correct leading behaviour for rj — > — oo. Of course, the 
particles with masses m n larger than twice the lightest mass m\ are unstable. The masses m n 
become dense above 2m as r\ — > — oo, and condensate to reproduce the continuum spectrum 
when h vanishes. The spectrum (|6.15|) was first obtained in [BS] through the study of the 
analytic structure in momentum space of the spin-spin correlation function for small magnetic 
field. Relativistic corrections to (|6.15]> have recently been obtained in |66] 135] . 

It is possible at this point to figure out the evolution of the particle spectrum of the Ising 
field theory in the plane of Figure 1. We have just described how the kinks living at rj = 
—oo are confined into a tower of bound states when the magnetic field is switched on. The 
number of stable particles (i.e. with m n < 2m\) decreases as r\ is increased from — oo to finite 
negative values. We expect three such particles to be left for small values of r\ (i.e. close to 
the magnetic axis), although five more particles above threshold happen to be stable at rj = 
due to integrability. The number of stable particles continues to decrease as rj increases, until a 
single particle is left at rj = +oo (the positive thermal axis). This scenario, which is an updated 
version of that originally proposed in [HSj, is supported by the numerical investigation of the 
spectrum of the Ising field theory [BUI IBB] . 

The correction to the mass of the particle when the magnetic field is switched on at r > 
can again be computed in perturbation theory. This time, however, the first order correction 
vanishes for symmetry reasons (the perturbing operator a couples only to an odd number of 
particles at r\ = +oo), and one has 



(the integral is evaluated at rj = +oo). In principle, one could think to evaluate the matrix 
element through a decomposition over intermediate particle states. Very recently, however, 
Fonseca and A. Zamolodchikov computed it generalising the techniques used for the spin-spin 
correlation function at h = 0. In this way they found |35j 




(6.16) 



where m is the mass at rj = +oo, and 





b = 21.52398.. 




within the normalisation (|5.2()]) of the operators. Since 



5£. 



vac — 



y J <Px(a(x)a(Q))c + 0(h*) 




7 / A h 2 + 0(h 4 ) 



(6.19) 
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we have 

|^ = -^(2 7 r)- 7 / 4 + 0(r / - 15 / 4 ) = -(11.66467..) +0(?/" 15/4 ), r? ^ +00 . (6.20) 

We mention, in conclusion, that |66j contains also a detailed study, which combines analytic 
expansions and numerical methods, of the analytic properties of the free energy of the Ising field 
theory as a function of temperature and magnetic field. 
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